Output of a pulsed atom laser 
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We study the output properties of a pulsed atom laser con- 
sisting of an interacting Bose-Einstein condensate (BEC) in 
a magnetic trap and an additional rf field transferring atoms 
to an untrapped Zeeman sublevel. For weak output coupling 
we calculate the dynamics of the decaying condensate popu- 
lation, of its chemical potential and the velocity of the output 
atoms analytically. 



The experimental breakthrough to Bose-Einstein con- 
densation with small numbers of atoms in magnetic traps 
has raised much interest in the properties of meso- 
scopic quantum gases. Bose-Einstein condensates with 
atoms in a single magnetic sublevel have been studied 
experimentally and theoretically. The recent experimen- 
tal and theoretical investigations of interference between 
two independent Bose-Einstein condensates convincingly 
proved their macroscopic coherence Moreover, the 

laser-like coherence of the atoms is preserved in the pres- 
ence of a matter-wave splitter based on rf-transitions 
pumping the atoms into untrapped magnetic sublevels 
[Q. These states are either strong-field seeking or have 
no magnetic moment at all, and leave the trap. Alter- 
natively, optical Raman transitions can be used for the 
transfer |^,^. Such schemes provide controllable output 
couplers for coherent atom lasers. 

In analogy to a laser one can distinguish between a 
cw laser, based on continuous refilling of the condensate, 
and a pulsed atom laser where the condensate is periodi- 
cally refilled and slowly released, similar to [Q . Whereas 
a continuous wave atom laser has been studied only theo- 
retically |^-|^ , current Bose-Einstein condensation exper- 
iments are limited to the pulsed mode of operation. The 
closest approximation of a cw atom laser by a pulsed one 
can be reached in the limit of a weak coupling rf field. In 
this case we are able to describe the decay of the trapped 
condensate and its energy width analytically. Previous 
calculations addressed the opposite limit of strong cou- 
pling by numerical calculations |^ or neglected the im- 
portant influence of atom-atom interactions [|o| . 

The output coupler consists of a monochromatic res- 
onant rf field of frequency Wrf transferring ^'^Na atoms 
in the F = 1 hyperfine state from the trapped m — —1 
into the untrapped m = and the repelled m — 1 mag- 
netic sublevels. For simplicity an isotropic harmonic trap 



potential V-i{r) = VoR + Afw|r72, V+i{v) = -V-i(r) 
and Vb(r) = are assumed while effects of gravity are 
neglected. 

The three coupled coherent matter waves are described 
by a three-component Gross-Pitaevskii equation (GPE) 
with resonant excitation in rotating wave approximation 
first studied for a generic two-level system in Ref. . 

In the following we adopt the point of view of sponta- 
neously broken gauge symmetry for a Bose gas initially 
at zero temperature. The system of equations for the 
macroscopic wave funtion ipmit) = e~^"^'^'^* {4>m{t)) in ro- 
tating wave approximation for m, m' € { — 1, 0, +1} now 
reads 
d ■- 
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Here, we have replaced the atomic density by the modu- 
lus of the wave function 



(n(r,i)) = ||V'(r,t)||^ 



(2) 



The coupling constant hi} — giJ,Bohi-\B\/ \/2 refers to the 
Rabi frequency due to the rf field. At zero magnetic 
field the symmetrized s-wave scattering matrix elements 
Umm' — ^T^T^^amm' for an elastic collision of a pair of 
atoms in the sublevels m,m' e { — 1,0, -1-1} are all nearly 
equal to a = 53ao, (oq is the Bohr radius), according to 
preliminary calculations of E. Tiesinga and P. S. Julienne 
pT[ . Since the steady state operation depends mainly on 
the initial condensate mean field, we assume in the fol- 
lowing a diagonal scattering matrix amm' = Smm'O, for 
simplicity. Consequently, the Hartree mean field poten- 
tial for each spin component is equal to the total atom 
density (n(r,t)) multiplied hy U — A-Kfi^a/M . 

The initial condition is chosen as the solution of the 
stationary GPE for the trapped condensate in the ab- 
sence of the rf field, i.e = in Eq. (|l|). In the Thomas- 
Fermi approximation it reads 



l^-i(r,0)|^ 



li + VoS-V-i{v) 
U 



(3) 



For a small coupling strength {Q, <C wt) the process of 
atoms leaking out of the resonance points is faster than 
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the Rabi oscillations. Therefore the coupling into state 
m = +1 can be neglected, since it is proportional to fl'^. 
In the following only the states to = — 1 and m — are 
considered. After switching on the coupling due to the 
rf field, initial oscillations die out quickly, because the 
untrapped atoms leak out of the trap within less than 
one Rabi cycle. Other condensate atoms move into the 
resonance area replacing the leaving ones. Eventually a 
quasi-stationary state is reached, i.e. the m = — 1 con- 
densate wave function decays slowly without oscillations 
while the atoms coupled out of the condensate are ex- 
pelled due to the mean-field potential and form a steady 
current. A numerical solution of the two-component 
GPE shows the uniform decay of the trapped condensate 
(cf . Fig. |l|) . The sum of external and mean field poten- 
tials is in a good approximation spatially independent 
inside the condensate so that a description of the decay 
in terms of the chemical potential ^{t) is appropriate. In 
the following we present an analytical calculation of the 
time dependent output intensity in this quasi-stationary 
regime, assuming a three dimensional isotropic harmonic 
trap. 

We first calculate the rate F of transitions from the 
condensate into the output. In the spirit of the Thomas- 
Fermi approximation we neglect the kinetic energy in a 
two-component GPE corresponding to Eq. (|l|) and solve 
for the output density distribution 



\Mr,t)\^ = 



4f]2 sin2[i^A2(r) +4172 , 



A2(r) +4172 



l^-i(r,0)|^ 
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The maximum amplitude of the Rabi oscillations is lo- 
cated at Trcs determined by the resonance condition 



nA(r,es) = 



(5) 



where hA{r) = hu-^f — V^i{r). Thus, the main contri- 
bution to the output coupling stems from a small shell 
around that resonance radius. 

The time-derivative of the density gives the density 
transition rate. The total transition rate is obtained by 
integrating the position dependent rate over the conden- 
sate volume. The transition rate is negligible outside a 
minute resonance shell and strongly peaked within that 
shell. Expanding the position dependent detuning A(r) 
to first order one obtains the time dependent transition 
probability 



T[t) 



1 



|^-i(rrcs,0)|2yo dt 
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, r sin[(7A'(r,cs)V+4r!2 1] 
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Jo(217i). 



(6) 



The calculation of the above rate V{t) does not account 
for the losses due to the leaving atoms. However, a rate 
equation allowing for these losses can be derived in the 
limit of weak coupling by using the perturbational rate 
{t < 7r/17) 



(Mtj2)3/2 



(7) 



with A'(r-rcs) = 2A(0)/r„s. 

Due to the spatial localization of the output coupling, 
the decay of the condensate population 



(8) 



N{t)= / d^r|V-i(r,i)|2 , N{0) = N^ 



depends solely on the density of the atoms around the 
resonance shell with radius r^s 



dN{t) 
dt 



~F|Vi-l(n.es,i)|' 



(9) 



In the quasi-stationary regime we assume the shape of 
the condensate density being equal to the Thomas-Fermi 
solution of the stationary GPE (cf. Fig. |l|) with a slowly 
varying atom number N(t), 



in 



N{t) = __[2/i(i)]5/2/(Mtj|)3/2. 
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(10) 



The condensate density at the resonance points is then 
given by 



\^Mrrcs,t)\^ = Ht)-nA{0)]/U. 



(11) 



Inserting this into the decay law (^ we obtain a non- 
linear differential equation for the decay of the chemical 
potential 



- nA{o) 



0, 



(12) 



where a = 3r(A/ti;^)'^/22 7/2^ ^ . Integration yields 



-2[;iA(0)]3/2artanhv//i/?iA(0) 



= at 



(13) 



With Eq. (|T^) this yields additionally the time evolu- 
tion for the number of trapped atoms N{t) as well as the 
flux and the velocity of the untrapped atoms. 

During the depopulation of the condensate the chemi- 
cal potential and the spatial extension of the condensate 
decrease until the resonance points lie on the surface of 
the shrinked condensate. At this point the flux out of the 
condensate vanishes, because 1-0-1 (?'rcs, i)P = 0, and the 
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chemical potential and the number of atoms remaining 
in the trap become constant in time 



^(oo) = ;iA(0). 



(14) 



As an example we chose a small value of A(0) and a 
trap frequency uj^ — 2t: x 106 Hz, the geometric mean 
of the values given in Ref. The resulting time evo- 
lution of the system variables according to the nonlinear 
differential equation is shown in Figs, ^(a) and (b). The 
chemical potential fj,{t) and the number of trapped atoms 
reach their steady state after roughly 25 seconds. Cor- 
respondingly the flux and the velocity of the untrapped 
atoms decrease to zero (cf. Figs. |2](c) and (d)). The 
corresponding calculations were also carried out in one 
dimension in order to compare the analytical expressions 
with numerical simulations of the full coupled GPE's. 
The results showed excellent agreement (less than 5 % 
deviation). 

The finite duration of the atom pulse ejected from the 
trap leads to a finite energy width of both the conden- 
sate and the output beam that is called here the natural 
energy width 5E of a pulsed atom laser, in analogy to the 
natural linewidth of spontaneously emitted photons. It 
can be calculated approximately by describing the initial 
stage of the output coupling process by an exponential 
decay of the condensate population 



^(O^-Fpop [N{t)~Ni^)]. 
The rate Fpop results in an energy uncertainty 



6E = hT, 



m(0)3/2 



(15) 



(16) 



Whereas the transition rate F does not depend on fi and 
Nq, the population decay rate Fpop does so due to its 
dependence on 'ijj-i{r^O). We thus find a dependence of 
the energy width on the condensate number 



SE 



N, 



3/5 



(17) 



The natural energy width becomes narrower for weaker 
coupling strength fl, i.e. a slower output coupling pro- 
cess. The same effect can be achieved by choosing a 
smaller detuning A(0) which causes the sphere of the 
resonance points to shrink towards the center of the con- 
densate. The energy width can further be reduced by 
starting out from a condensate with a large population 
No- 

The energy width SE can also be understood as a ve- 
locity width 6v — SE/ (Mv) of the untrapped atoms leav- 
ing the condensate. The parameters given in Fig. ^ result 
in (v) — 1.31cm/s for the velocity outside the condensate 
and Sv/ (v) w 10~^ for the relative velocity width. 



An additional energy width is imposed by the tempo- 
ral decay of the chemical potential nit). It leads to a 
decrease of the output velocity according to Mv{t)^/2 — 
fi(t) — 7iA(0) thereby implying a frequency chirp of the 
output beam. This frequency chirp can be compensated, 
however, by imposing a chirp on the frequency cjj-f of 
the rf field such that the output velocity rather than the 
detuning A(0) becomes constant. 

The output pulse can last for up to the order of 100 sec- 
onds such that causes of phase fiuctuation which would 
also influence a cw-atom laser have to be considered as 
well. Additional broadening of the energy width of the 
output beam is caused by thermal excitations of the con- 
densate wavefunction and technical noise of the output 
coupling mechanism such as fluctuations of the conflning 
magnetic fleld. Thus, the natural energy width calculated 
above has to be understood as a lower limit. A compre- 
hensive theoretical description of the fluctuations of the 
condensate wavefunction, particularily of its phase, is be- 
yond the scope of this Letter. 

It has been shown in |^ that the temporal decay of 
phase correlations transforms into a decay of the spatial 
coherence of the output beam along the mean classical 
trajectory that is being passed during the according cor- 
relation time of the condensate phase. 

At present, random variations of the bias magnetic 
field at the mGauss level are among the most important 
experimental limitations to the coherence properties of 
the output, prior to fluctuations of the condensate wave- 
function. They might contribute to a phase diffusion of 
the output beam in the 100 Hz range p2[ . 

Pulsed atom lasers will play an important role in cre- 
ating coherent matter waves as the size of the condensate 
may be largely increased in future experiments. The re- 
sults of this Letter are based on the solutions of the cou- 
pled Gross-Pitaevskii equations and allow us to extract 
the relevant properties of the output analytically. The 
velocity of the untrapped atoms depends on the slowly 
decaying population of the trapped condensate fraction, 
leading to a slow chirp of the output frequency. While 
this frequency chirp might be compensated by a variation 
of the rf-fleld frequency, other sources of spectral width 
like fluctuations of the confining magnetic field are surely 
more difficult to overcome. 
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FIG. 2. Time evolution in the Thomas-Fermi approxima- 
tion: (a) chemical potential, (b) number of particles in the 
trapped m = —1 state, (c) particle flux in the untrapped 

m — state, (d) velocity of atoms leaving the conden- 
sate. The parameters are A'^o = 5 x 10®, A(0) = 3500 s"\ 
n = 20s"^ 



FIG. 1. Uniform decay of the condensate density in terms 
of the effective potential V^-i(r) — VoB + L''|'^^-i(r)|^ for small 
field strength « = 12 s-\ iVo = 5 x 10^ and A(0) = 3100 s"^ 
in a one-dimensional situation. The lower density is reached 
after 670 ms. 
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